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SYMBOLS 

The  following  notation  is  used  in  the  ANALYSIS: 
speed  of  sound 


a 


solution  of  hypergeometric  equation  (27) 


K 


coefficient  in  Chaplygin’s  equation,  defined  hy 
equation  (20) 


K* 


step  function  approximating  function  K 


K* 


function  defined  "by  equation  ( 38) 


M 


Mach  number 


n 


normal  to  the  curve  T  pointing  toward  the  domain  A 


[cr 


speed 

critical  speed 


‘max 


maximum  value  of  q 


arc  length  measured  on  T 


length  of  f 


function  of  speed  (see  equation  (l8)) 


^m 


ax 


value  of  t  for  M  =  00 


Preceding  Pag^S  Blank 
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* 


* 


* 


V 


u,v 

x,y 

X(s),Y(s) 

7 

r 

rh 

A 

e 

p 

Per 

T 

Tcr 

0 

♦ 

CO 

A; 

0 


components  of  velocity  vector 
coordinates  in  physical  plane 
coordinates  of  points  on  T 
exponent  in  pressure— density  relation 
curve  in  physical  plane  along  which  M  =  1 
image  of  r  in  hodograph  plane 
domain  in  which  a  subsonic  flow  is  given 
angle  between  velocity  vector  and  x— axis 
density 

critical  density 

variable  defined  by  equation  (28) 
value  of  t  for  M  =  1 

velocity  potential 

stream  function 

angle  between  tangent  to  f  and  x— axis 

used  over  symbol,  value  on  T 
used  as  subscript,  value  for  q  =  0 


L 
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The  following  notation  is  used  in  the  appendix: 

end  points  of  segment  along  which  Cauchy  data  are  given 
A,B,T,N,T' ,N' ,T",N"  constants  defined  in  section  3  of  the  appendix 
D(a^-b)  domain  defined  in  section  1  of  the  appendix 

K  coefficient  in  equation  (Al) 

functions  approximating  function  K 
complex-valued  functions 


x,y 

Y(n)^(n) 


independent  variables 

functions  defined  by  equations  (Al8) 


z,(, 

Zn,i*Z^n^ 

AV 

S(K) 

T,V 

9 

♦ 


complex  variables 

"formal  powers"  defined  in  section  4  of  the  appendix 
constants  defined  by  equations  ( A3^-) 
matrix  defined  by  equation  (A9) 

given  initial  values  of  \Jr  and  its  normal  derivative 
function  connected  with  by  equation  (AT) 
dependent  variable  in  equation  (Al) 
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.  ANALYSIS 

1.  The  Continuation  Theorem 


The  main  result  of  this  investigation  reads  as  f ollovs : 

Given  a  domain  A  in  the  x,y^-plane,  whose  boundary  contains 
a  continuously  curved  arc  T  (see  fig.  l);  given  a  subsonic 
potential  gas  flow  in  A;  and  if 

(a)  the  Mach  number  of  the  flow  on  r  is  equal  to  unity, 

(b)  the  components  of  the  velocity  vector  and  their  first- 
order  derivatives  are  continuous  on  T, 

(c)  the  nonnal  derivative  of  the  speed  on  r  is  different 
from  zero,  and 

(a)  no  streamline  of  the  flov  cuts  r  at  a  right  angle j 
then  the  flow  can  be  continued  across  any  subarc  of  r  as  a 
potential  supersonic  flow  without  weak  discontinuities  and  this 
continuation  is  uniquely  determined. 

This  theorem  contains  no  general  statement  concerning  the  extent 
of  the  supersonic  region  obtained;  however,  in  each  special  case  the 
method  of  continuation  described  as  follows  determines  the  extent. 

By  a  compressible  potential  flow  without  weak  discontinuities  is 
meant  a  flow  described  by  two  twice  continuously  differentiable 
functions  0(x,y)  and  ^(x,y)  (potential  and  stream  function) 
satisfying  the  differential  equations 
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where  pQ  is  the  stagnation  density  end  p  the  density  given  "by 


Here  7  is  the  exponent  in  the  pressure-density  relation}  a^  is  the 
speed  of  sound  at  a  stagnation  point}  and  q,  the  speed,  is 


It  is  well  known  that  in  the  supersonic  region  flows  with  weak 
discontinuities  ( discontinuities  in  the  second  derivatives  of  the . 
velocity  potential)  are  possible*  In  fact,  there  is  no  good  physical 
reason  to  assume  that  in  a  supersonic  flow  the  velocity  components 


possess  partial  derivatives.  (The  continuity  equation  and  the  condition 
for  the  existence  of  a  potential  may  he  expressed  in  terms  of  integral 
relations.)  Under  precisely  which  conditions  flows  with  weak  disconti¬ 
nuities  can  he  obtained  hy  continuing  subsonic  flows  across  the  sonic 
line  seems  to  he  a  very  important  hut  rather  delicate  question}  it  is 
intimately  tied  up  with  the  question  of  the  necessity  of  conditions  (h) 
to  (d) . 

If  condition  (c)  is  violated  along  an  arc  of  T,  then  the  velocity 
vector  is  constant  along  this  arc.  This  follows  from  equation  (8) .  It 
seems  very  probable  that  in  such  a  case  the  velocity  vector  must  be 
constant  everywhere.  Whether  the  continuation  theorem  remains  true  if 
condition  (c)  is  violated  at  isolated  points  remains  an  open  question. 
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Condition  (d)  is  probably  essential  for  the  uniqueness  of  the 
continuation.  If  at  some  point  P  of  T  the  streamline  is  normal 
to  P,  then  the  Mach  lines  at  this  point  are  tangential  to  r,  for 
at  P  the  Mach  angle  is  90°.  Such  a  situation  arises  in  a  synmetrical 
Laval  nozzle  if  the  transition  from  subsonic  to  supersonic  speeds  occurs 
along  a  line  extending  across  the  nozzle  (flow  of  Mayer  type,  see 
fig.  2).  In  a  remarkable  paper  (reference  2)  Frankl  investigated  this 
case  and  showed  that  the  continuation  of  the  subsonic  flow  into  the 
supersonic  region  is  not  uniquely  determined.  He  showed  also  that  in 
this  case  the  supersonic  flow  may  exhibit  weak  discontinuities. 


* 


4 


2.  Reduction  to  a  Cauchy  Problem  in  the  Holograph  Plane 
In  order  to  prove  the  theorem,  let  the  equation  of  the  curve  F  be 


x  =  X(s) 
y  =  Y(s)  ► 
0  <  8  < 


(5) 


where  s  is  the  arc  length  of  F  measured  in  the  direction  of 
traversing  F  with  the  domain  A  to  the  left.  By  hypothesis  X"(s) 
and  Y"(s)  exist  and  are  continuous,  and  X*^  +  Y*^  =  1.  It  follows 
that  the  angle 


m(s)  =  tan-1  (6) 

X*(s) 


between  the  tangent  to  T  at  x  =  X(s),.  y  =  Y(s)  and  the  positive 
x-direction  is  a  continuously  differentiable  function  of  s.  (See 
fig.  1.) 


a 


u 


Set 
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The  value  of  any  function  F(x,y)  at  x  =  X(s),  y  =  Y(s)  will  be 
denoted  by  F(s).  Thus  0(s)  =  @[x(s) ,  Y(s)J.  Nikolskii  and  Taganov 
(reference  3)  showed  that 


0*( s)  =  COS2(e  -  CD) 

V 2ao2  ^ 


(8) 


where  b/'dn  denotes  differentiation  in  the  direction  of  the  normal 
to  r  pointing  toward  the  domain  A.  (This  formula  follows  easily 
from  the  equation  of  motion  written  with  q  and  0  as  independent 
variables.)  Since  M  =  1  on  T  and  M<1  in  A  (M  =  q/a  being 
the  Mach  number)  it  follows  that  dq /bn  <  0,  so  that  by  hypotheses  (c) 
and  (d) 


0'( s)  <  0 


(9) 


It  is  no  loss  of  generality  to  assume  that 


2*  >  0(0)  >  0(S)  >  0  (10) 


If  this  condition  were  not  satisfied,  T  could  be  divided  into  a 
finite  number  of  overlapping  arcs  along  each  of  which  0  would  change 
by  less  than  2it.  Inequalities  (9)  and  (10)  imply  that  F  has  a  one— 
to-one  image  in  the  holograph  plane,  which  is  the  circular  arc  F^ 
given  by 


<1  =  <lcr,  0(S)^0<0(O)  (11) 


where  qcr,  the  critical  speed,  is 


y 
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■* 


By  the  definition  of  the  potential  and  stream  function, 


=  —  (u  di  +  y  dy) 
a0 


(12) 


dt  =  (— v  dx  +  u  dy) 


p0a0‘ 


On  r. 


* 


~  =  qcr  cos  & 

v  =  qcr  sin  6 


dx  =  ds  cos  cd 

dy  =  ds  sin  (D 


so  that 


* 


^(s) 


■  (7Tt)1/2j[' 


cos  (0  —  cd)  ds 


+  (s)  =  ~(y  +  J  sin  (§”  -00)  ds 


(13) 


/\4  'V 

These  equations  show  that  ^(s)  :^and  ^(s)  are  twice  continuously 
differentiable .  By  hypothesis  ?(s)  is  a  continuously  differentiable 
function.  By  virtue  of  inequality  (9)  the  inverse  function 


% 
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3  =  s*(e) 

exists  and  hence  is  continuously  differentiable.  It  follows  that 


0*(e)  =  0 

\ 

js*(9)] 

► 

\jf*(  9)  =  \jr 1 

s*(9)] 

(iM 


are  continuously  differentiable  functions  of  0,  9(s)  <  9  <  9(0) . 


Assume  now  that  the  flow  in  A  has  been  continued  across  T . 

Since 


d(u,v)  _  qd(q,9) 
8(x,y)  8(x,y) 


it  follows  that  on  T 


8(u,v) 

8(x,y) 


cos^(9  —  cd) 


(15) 


In  fact,  on  T 

8(q,9)  8q  89  8q  89 
8(x,y)  8s  8n  8n  8s 

and  8q/8s  =  0,  whereas  89 /8s  is  given  by  equation  (8).  By 
hypotheses  (c)  and  (d)  expression  (15)  cannot  vanish.  Hence 


Mgall  <  0 

8(1,  y) 


for  all  points  sufficiently  close  to  F.  It  follows  that  a  sufficiently 
small  neighborhood  of  any  subarc  of  X  has  a  one-to-one  image  in  the 
holograph  plane,  that  is,  a  neighborhood  U  of  a  subarc  of  In  the 
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holograph  plane  the  functions  0  and  i  satisfy  the  well-known 
Chaplygin  equations 

\ 

50  ^PQ  5i 
50  P  5q 

?  (16) 


M  1  -  M2  PO  5t 
5q  4  p  80^ 


while  on  they  satisfy  the  initial  conditions 

0(qcr,0)  =  0*(0)^ 

► 


t(qcr*0)  =  **(0)  J 


(17) 


for  0(S)  <  0  <  0(0). 

3.  Existence  and  Uniqueness  of  the  Solution 
in  the  Hodograph  Plane 

In  the  preceding  section  the  continuation  of  a  subsonic  flow 
across  the  sonic  line  has  been  reduced  to  the  initial— value  problem, 
equations  (l 6)  and  (17)«  Introduce  the  new  independent  variable 

t  =  p  £-k  (18) 

iKcr  ^ 

(Note  that  q  =  0  corresponds  to  t  =  — <»,  q  =  qcr  to  t  =  0, 
and  q  =  qmax  =  Sq  'J& /{y  —  l)  to  some  finite  positive  value  of  t.) 

In  the  0,1>-plane,  0  and  i]r  satisfy  the  equations 
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where 


M 

So  St 

> 


M,^(t  fit 

St  So  J 


(19) 


K(t)  =  L=J£  (20) 

2 

P 


The  initial-value  problem,  equations  (l6)  and  (17),  goes  over  into 
a  Cauchy  problem  for  equations  (19).  The  theory  of  this  Cauchy  problem 
is  given  in  the  appendix.  From  the  theorems  stated  in  section  3  of  the 
appendix  it  follows  that  the  initial— value  problem,  equations  (l6) 
and  (l7) y  possesses  a  unique  solution  in  the  domain  bounded  by  the 
.arc  rh  and  two  characteristics  of  equations  (l6) .  As  is  well  known, 

these  characteristics  are  epicycloids.  In  this  domain  (see  fig.  3)> 
which  is  called  the  "characteristic  triangle  determined  by 

the  functions  0  and  \[f  will  be  continuously  differentiable,  provided 
that  0*(6)  and  ty*(e)  have  this  property. 

Remark.—  Note  that  if  ^  (considered  as  a  function  of  Q  and  .t) 
is  twice  continuously  differentiable,  then  equations  (19)  iinply  that 
satisfies  Chaplygin fs  equation 


de2  St2 


(21) 


* 


* 


In  this  case  the  initial— value  problem  considered  may  be  formulated  as 
a  Cauchy  problem  for  equation  (21) :  ^ 


ilr(0,O)  =  t(0) 

/  ' 

9 

tt(e,o)  =  v(e) 
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where 


T(0)  =  t*(e) 
v(e)  =  0**(e) 


In  general,  however,  *vjr  need  not  be  a  twice  continuously 
differentiable  function  in  the  supersonic  part  of  the  hodograph  plane. 

4.  Existence  and  Uniqueness  of  the  Solution 
in  the  Physical  Plane 

In  order  to  complete  the  proof  of  the  continuation  theorem,  it 
must  be  shown  that  the  functions  0  and  i]r  obtained  in  the  preceding 
section  can  be  transferred  (in  a  unique  way)  to  the  physical  plane 
(more  precisely,  to  a  part  of  the  physical  plane  adjacent  to  r  and 
lying  outside  A)  and  that  in  the  x,y^-plane  these  functions  possess 
continuous  derivatives  of  the  second  order  which  assume  the  "correct" 
values  along  I\ 

The  functions  0  and  obtained  from  the  Cauchy  problem  coincide 
along  F  with  the  potential  and  stream  function  of  the  flow  given 
in  A  and  transferred  to  the  hodograph  plane  by  means  of  the  mapping 


(22) 


It  follows  from  equations  (12)  that  the  inverse  mapping  is  given  by 


x  +  iy  =  Q.Q  J  fii^*|d0  +  i-2.  dty 


u  =  u(x,y) 

v  =  v(z,y) 


(23) 


Ik 
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or. 


Since  the  functions  0  and  \]r  have  "been  continued  into  the 
supersonic  part  of  the  hodograph  plane,  this  integral  is  meaningful 
even  there.  Moreover,  it  is  path— independent.  Under  the  assumption 
that  the  partial  derivatives  0ee,  00q,  ...  exist  and  are  continuous, 

this  follows  immediately  from  Chaplygin's  equations  (16).  But,  in  the 
case  considered,  this  assumption  cannot  he  made.  Consider,  however,  the 
mapping  from  the  hodograph  plane  to  the  0, ty— plane.  By  virtue  of 
equations  (l6) 


=  /Mf  +  £1  /Mf 

3(0,q)  p0q  W  p0(l  -M2)  \3q/ 


It  follovs  (for  instance,  from  relations  (9)  and-  (13))  that  this 
Jacobian  is  different  from  zero  in  some  neighborhood  of  the  arc 

Hence  the  mapping  from  the  hodograph  to  the  0,^— plane  is  locally  one- 
to-one,  so  that  it  is  sufficient  to  establish  the  path-independence  of 
the  integral,  equation  (23)>  in  the  0,^— plane •  This  leads  at  once 
to  the  conditions 

8  /cos  9\  b_  (p o  sln  eV 
3t\  q  )  80  V  pq  / 

8  /sin  8_  ( PQ  cos  0 
8t\  q  /  30  \  pq 


Using  the  relation 
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these  conditions  may  he  transformed  into 

"X 

80  _  8q 

^  P0<1  3* 

> 

Be  PpC1  -  m2)  Bg 
St  pq.  80 


(24) 


Interchanging  the  dependent  and  independent  variables ,  it  is  seen  that 
system  (24)  is  equivalent  to  equations  (l6).  Thus  equations  (l6)  insure 
the  path-independence  of  the  integral,  equation  (23),  even  without  any 
assumption  about  the  existence  of  second— order  derivatives  of  0 
and  t* 

Thus  equation  (23)  defines  a  mapping  also  in  the  supersonic  part 
of  the  hodograph  plane.  The  Jacobian  of  this  mapping  is  easily 
computed  to  be  equal  to 


— 2!  It!  _  m2)  2°i(2+f  *  (M) 

u,v)  q^  P  [_  p2  \Be/  \B QV 


(25) 


so  that  for  q  =  qcr 


7-2 


8(x,y) 

8(u,v) 


f i)"1  ^ 


But  by  equations  (l^) 


a*' (9)  =  M  Ml  =  (  2  \1//g  cos(e  -  co)  , 
r  v  ds  d9  \7  +  lj  ^'(s)  f 


so  that 


$(* ill  4  0 

8(u,v)  f 
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on  the  arc  defined  by  equation  (ll)  and  hence  also  in  some  neighborhood 
of  any  subarc  of  this  arc.  It  follows,  that  such  a  neighborhood  is 
mapped  onto  a  neighborhood  of  a  subarc  of  F  in  a  one— to— one  manner. 

The  functions  0  and  i  obtained  by  solving  the  initial— value  problem, 
equations  (l6)  and  (17),  can  now  be  transferred  to  the  x,  y^-plane  by 
means  of  the  mapping  defined  by  equation  (23). 

The  derivatives  of  0  and  i|r  in  the  x,  y-plane  can  be  computed 
by  using  equations  (l6)  and  (23).  (See  reference  4.)  This  direct 
though  lengthy  computation  yields  exactly  equations  (12) .  Thus  0 
and  i|r  satisfy  equations  (l).  Finally,  equation  (23)  shows  that  x 
and  y  are  continuously  differentiable  functions  of  u  and  v,  so 
that  u  and  v  are  continuously  differentiable  functions  of  x 
and  y.  Thus  0  and  i|f  possess  continuous  partial  derivatives  of  the 
second  order  in  the  x, y-plane.  The  proof  is  now  complete. 


5.  Computation  of  the  Solution 

The  theorems  proved  in  the  appendix  implicitly  contain  methods  for 
the  effective  computation  of  the  solution.  Thus  theorem  2  shows  that 
the  solution  of  the  Cauchy  problem  may  be  represented  as  in  infinite 
series  of  particular  solutions  of  the  Chaplygin  type.  (The  idea  of  the 
following  proof  is  already  contained  in  a  remark  by  Frankl  concerning 
the  Tricomi  problem;  see  reference  5 • ) 

Chaplygin  (reference  6)  showed  that  a  particular  solution  of  the 
second— order  equation  obtained  from  equations  (l6)  by  eliminating  0 
can  be  found  by  setting  (for  any  positive  constant  a) 


+  =  eia®(q/ (3jnax)a®'a  (^Amax2)  ^ 

where  q„  '  is  the  maximum  speed  possible  for  the  given  gas: 


and  F<x(t)  is  any  solution  of  the  hypergeometric  equation 
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T(1  -  T)Fa"(T)  +  [l  +  a  -  (aa  +  ba  +  !>]  V(t )  -  aabaFa(T)  =  0  (27 ) 

the  constants  and  ha  "being  determined  by  the  equations 


ao^a  " 


a(a  +  l) 
2(7  - 1) 


It  is  seen  that  the  potential  corresponding  to  equation  (26)  is 
given  by 

<j6  =  -ielafi(p0/p)(q/qmax)a  [Fa(  t)  _  (2T/a)Fa*  (  t)J 

where 

T  =  (28) 

Let  Fqj(t)  and  F(Xh(t)  denote  the  particular  integrals  of 
equation  (27)  which  for 


Ijnax 


satisfy  the  initial  conditions 
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0aj  =  -  0  for  <1  =  9cr^  -5  =  1 


"  °*  "  1  for  4  =  <lcr>  >5  =  II 


(31) 


Consider  nov  the  Cauchy  problem,  equations  (l6)  and  (17),  assuming 
for  the  sake  of  simplicity  that 


e(s)  =  -t 

0(0)  =  T 
0  <  T  <  n 


Set 


n  -  a* 
an  -  T 


(32) 


and  expand  the  initial  data  0*(0)  and  ty*(0)  in  Fourier  series 


1*(0)  =  /  an  cos  ctn©  +  bn  sin  0^0 
n=0 


\|r*(0)  =Xlcn  003  “n9  +  In  sIn 

n=0 


(33) 


T  <  0  <  T 


From  the  way  the  particular  solutions,  equations  (30),  have  "been 
constructed  it  follows  that  the  series 


1 


20 


NACA  TN  2058 


00 

0  =  a0  +1>~  (a^  cos  a^e  +  bn  sin  On^/AonlM  + 
n=l 


.  « 


*>~  (cn  sin  otn©  -  dn  cos  a^A^ntr) 
n=l 


> 


(cn  cos  On©  +  dn  sin  an©)BaniI(T) 


m 


n=l 


(-an  sin  0^0  +  t>n  cos  “n^Bot^I^) 


y 


represent  a  formal  solution  of  the  initial— value  problem  considered. 


Theorem  2  of  the  appendix  not  only  shows  that  the  series  converge 
absolutely  and  uniformly  (within  the  "characteristic  triangle")  to  the 
desired  solution  but  also  permits  an  estimate  of  the  error  committed  by 
replacing  the  infinite  series  by  finite  partial  sums.  Theorems  4  and  5 
of  the  appendix  contain  similar  statements  concerning  the  convergence 
of  the  differentiated  series. 

Remark  1.-  While  the  theorems  in  the  appendix  refer  to  an 
equation  of  the  form  of  equation  (2l),  it  is  clear  that  these  theorems 
imply  similar  ones  for  system  (l6).  It  may  be  noted  that  the  convergence 
of  series  (34)  was  proved  by  a  general  method;  no  special  properties  of 
the  hypergeometric  functions  were  used. 


Remark_2.—  It  should  be  noted  that  if  the  arc  Th  on  which  the 
initial  values  are  given  exceeds 


* 


then  the  epicycloids  drawn  from  the  extremities  of  the  arc  do  not 
intersect.  In  this  case  the  characteristic  triangle  must  he  replaced 
by  the  domain  bounded  by  the  initial  arc/  the  two  epicycloids,  and  an 
arc  of  the  circle  q  =  qmax.  (See  fig.  4.) 
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6.  Computation  of  the  Solution. hy  a  Second  Method 

The  expansion  in  a  Chaplygin  series  may  not  be  the  most  efficient 
way  of  computing  the  solution  of  the  initial— value  problem  in  the 
hodograph  plane,  since  it  requires  the  preliminary  computation  of  the 
functions  A&  and  Ba,  and  because  of  the  possibility  of  slow 
convergence.  Another  method  is  suggested  by  theorems  3  and  6  of  the 
appendix. 

Let  the  equation  satisfied  by  the  stream  function  be  transformed 
into  the  form  of  equation  (23).  Divide  the  interval  0  <  t  <  t^^ 
into  subintervals 


and  set 


t  y <  t  <  ty 

V  =  1,  2  .  .  .  n 


(35) 


0  =  to  <  ti  <  .  .  .<  tn  -  tmaT 


K*(t)  =  min  K(t' ) 


(36) 


tv_2_  <  t  <  tv  ,  tv_l  ^  t  *  <  ty 

The  equation 


K*(t)*e0  +  1rtt  "  0  (3T) 


is  an  approximation  to  equation  (21)  .  Since  K*(t)  is  a  piecewise 
constant  function  ( step  function) ,  it  is  equivalent  to  the  one¬ 
dimensional  wave  equation  in  each  strip  tv  <  t  <  ty+-j_.  A  solution  of 
equation  (37)  will  be  required  to  possess  continuous  partial  derivatives 
of  the  first  order.  The  potential  0  corresponding  to  a  given 
solution  of  equation  (37)  “will  be  continuous  but  its  first  derivatives 
will  possess  discontinuities.  The  characteristics  of •  equation  (37)  are 
given  by  the  equations 
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Q 


-/4 


dt  =  Constant 


+  f  >/-£*(  tj 


dt  =  Constant 


vhich  represent  polygonal  arcs  in  the  6 , tr-plane .  (See  figs.  5  and  6.) 

For  equation  (37)  the  solution  of  the  Cauchy  problem, 
equations  (24),  is  a  trivial  matter.  This  solution  may  be  constructed 
by  repeated  use  of  the  classical  D'Alembert  foimula  (as  described  in 
detail  in  section  6  of  the  appendix) .,  Another  possibility  is  to.  apply 
the  method  of  the  preceding  section. 

Theorems  3  and  6  assert  that  for  a  sufficiently  fine  subdivision, 
equations  (35),  the  solution  of  equation  (37)  will  be  an  arbitrarily 
close  approximation  to  the  solution  of  equations  (24)  with  the  same 
initial  values. 

Instead  of  using  the  approximating  step  function,  equation  (36), 
one  may  use  the  function 

K*(t)  =  max  K(t»)  (3? 


tv-1  <  t  <  tv ,  tv_x  £  t-‘  £  tv 

In  this  case  K*(t)  =  0  for  0  <  t  <  tx,  so  that  in  this  strip  +  is 
a  linear  function  of  t. 

In  some  cases  it  might  he  preferable  to  use  as  the  approximating 
coefficient  a  piecewise  linear  function.  In  this  case  the 
approximating  equation  will  be  of  the  type  of  the  Darboux— 1 Tricomi 
equation 

^xx  +  +yy  =  °  •  (39) 

in  each  strip  tv  <  t  <  tv+1.  The  solution  of  the  Cauchy  problem  for 

the  Darboux— Tricomi  equation  can  be  expressed  either  by  integral 
formulas  (see  reference  7)  or  by  infinite  series  involving  Bessel 
functions . 
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7.  Application  to  Transonic  Flovs  past  Obstacles 

The  theorem  stated  in  section  1  is  of  a  purely  local  character. 
While  in  the  hodograph  plane  the  flow  can  be  continued  throughout  the 
whole  characteristic  triangle,  the  transition  to  the  physical  plane  may 
be  possible  only  in  the  immediate  neighborhood  of  the  sonic  line  T . 

This  may  be  due  either  to  the  existence  of  limiting  lines  or  simply  to 
the  fact  that  a  simply  covered  domain  in  the  hodograph  plane  corre¬ 
sponds  to  a  multiply  covered  domain  in  the  physical  plane. 

In  some  important  cases,  however,  it  can  be  asserted  that  the 
continuation  of  the  subsonic  flow  by  the  method  described  in  sections  2 
to  5  yields  the  whole  supersonic  region  which  is  of  relevance  to  the 
problem  considered. 

Consider  a  transonic  flow  past  a  closed  body  with  a  subsonic  stream 
Mach  number  (and  without  shock  waves).  The  supersonic  regions  Sp 
and  S2  (see  fig.  7)  are  bounded  by  solid  walls  and  transition  lines  • 
along  which  M  =  1.  For  a  supersonic  region  thus  bounded,  Nikolskii 
and  Taganov  (reference  3)  proved  an  important  theorem  stating  that  the 
region  S  possesses  a  one— to— one  image  in  the  hodograph  plane. 

In  this  proof,  Nikolskii  and  Taganov  make  use  of  the  remark  that 
nil  Mach  lines  in  S  must  possess  points  in  common  with  the  transition 
line  (the  line  M  =  l) .  While  they  give  no  formal  proof  of  this  fact, 
a  proof  can  be  supplied  without  difficulty.  The  domain  S  is  simply 
covered  by  a  family  of  smooth  curves,  the  streamlines  of  the  flow,  and 
the  solid  wall  bounding  S  belongs  to  this  family.  A  Mach  line  having 
no  points  in  common  with  the  transition  line  would  have  to  originate 
and  terminate  at  the  solid  wall.  A  simple  argument,  essentially  equiva¬ 
lent  to  Rolle's  theorem,  shows  that  such  a  Mach  line  would  be  tangent 
to  a  streamline  at  at  least  one  point  P.  But  this  is  impossible 
since  at  P  the  Mach  angle  would  have  to  be  0,  which  corresponds 
to  M  =  00,  p  =  0. 

Consider  now  the  hodograph  image  of  S.  It  is  bounded  by  the 

arc  rh  of  the  circle  q  =  qcr  (corresponding  to  the  transition  line 
in  the  physical  plane)  and  the  curve  Wh  (the  hodograph  image  of  the 
solid  boundary) .  Assume  that  Sh  is  not  contained  within  the  character 
istic  triangle  determined  by  which  is  bounded  by  and  two 

characteristics  of  the  hodograph  equations  (i.e.,  two  Busemann 
epicycloids) .  Then  at  least  one  of  these  Busemann  epicycloids 
intersects  Wh  at  two  points,  and  Sh  contains  an  arc  of  an 
epicycloid  which  has  no  points  in  common  with  This,  however,  is 
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impossible  since  the  epicycloids  are  the  hodograph  images  of  the  Mach 
lines.  Thus  is  contained  within  the  characteristic  triangle.  By 

virtue  of  the  previous  results,  this  implies  the  following  uniqueness 
theorem: 


Two  transonic  (subsonic  at  infinity)  flows  past  an  obstacle 

are  identical  if  they  are  identical  in  the  subsonic  region. 

Similar  theorems  are  true  for  flows  past  curved  walls  and  for 
Taylor  type  flows  in  nozzles  (figs.  8  and  9)*  In  the  case  of  a  Mayer 
type  nozzle  flow,  however,  the  subsonic  flow  does  not  determine  the 
supersonic  flow  .uniquely. 

The  preceding  remarks  are  of  interest  in  connection  with  various 
attempts  to  construct  transonic  flows  by  the  " correspondence  method." 

The  correspondence  method  was  first  used  by  Chaplygin  for  jet 
problems  (where,  however,  it  is  identical  with  the  solution  of  the  so- 
called  direct  problem) .  The  correspondence  method  may  be  expressed  in 
various  analytical  forms;  however,  except  for  the  very  original  method 
of  Bergman  (reference  8)  the  basic  procedure  is  always  the  same,  as  is 
clearly  pointed  out  by  Gelbart  (reference  9)-  It  consists  of  associ¬ 
ating  with  a  solution  of  the  Cauchy-Riemann  equations  (representing  an 
incompressible  flow)  a  solution  of  Chaplygin* s  hodograph  equations. 

In  some  cases  the  resulting  compressible  flow  is  of  the  same  general 
character  as  the  initial  incompressible  flow.  To  date,  the  most 
detailed  application  of  this  method  to  flows  past  obstacles  is  that 
given  by  Tsien  and  Kuo  (reference  10).  The  continuation  of  a  sub¬ 
sonic  flow  across  the  sonic  line  invariably  occurs  in  the  application 
of  the  correspondence  methods  to  transonic  flows. 

While  the  mathematical  problems  connected  with  the  correspondence 
method  are  as  yet  largely  unsolved,  the  results  of  the  present  investi¬ 
gation  show  that  once  the  solution  of  Chaplygin1 s  equation  in  the 
subsonic  part  of  the  hodograph  plane  is  determined,  the  continuation 
of  the  solution  in  the  supersonic  region  is  uniquely  determined  and  can 
be  effectively  computed.  If  this  continuation  does  no't  yield  the 
desired  result  in  the  physical  plane,  the  choice  of  the  solution  in  the 
subsonic  domain  must  be  abandoned. 


CONCLUDING  REMARKS 


The  continuation  of  a  given  subsonic  flow  into  the  supersonic 
region  has  been  discussed.  The  corresponding  problem  in  the  subsonic 
region  is  of  little  interest,  since  in  the  subsonic  region  the 
differential  equations  are  of  elliptic  type  with  analytic  coefficients 
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and  have  analytic  solutions.  Thus  if  the  velocity  potential  0  is 
known  up  to  a  line  T,  the  continuation  of  0  across  T  is  a 
problem  in  analytic  continuation.  It  seems  imp r oh able  that  any 
interesting  and  general  statements  can  he  made  concerning  this  case. 

In  the  supersonic  region  the  equations  are  of  hyperbolic  type  and  the 
continuation  of  a  given  flow  across  a  line  r  reduces  to  the  standard 
Cauchy  problem.  The  theory  of  this. problem  is  known  (see  reference  11, 
p.  326,  and  the  literature  quoted  therein).  It  would  be  very 
interesting,  however,  to  obtain  results  concerning  a  problem  converse  to 
the  one  considered  herej  that  is,  assume  that  a  supersonic  flow  is 
known  up  to  the  line  T  on  which  the  Mach  number  is  equal  to  1,  and 
determine  under  which  conditions  the  flow  can  be  continued  into  the 
subsonic  (elliptic)  region. 


Syracuse  University 

Syracuse,  N.  Y.,  September  1,  19^7 
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APPENDIX 

THE  CAUCHY  PROBLEM  FOR  CHAPLYGIN'S  EQUATION 


This  appendix  contains  the  proofs  of  the  theorems  used  to 
establish  the  preceding  results .  For  the  convenience  of  the  reader,  the 
notation  standard  in  the  theory  of  partial  differential  equations  is 
used.  Note  that  x  and  y  no  longer  refer  to  the  physical  plane. 
Equation  (Al)  is  essentially  the  same  as  equation  (2l) . 

1 .  Introduction 

The  aim  of  this  investigation  is  to  integrate  the  differential 
equation 


+  *„  -  0 


where 

K(0)  =  0" 
K(y)  <  0  ► 
for  y  <  0 

y 

under  the  initial  conditions 


(Al) 


(A2) 


t(x,0)  =  t(x) 


ty(x,0)  =  v(x) y 
a  i  xl  b 


(A3) 
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t  and  V  toeing  given  functions.  The  solution  is  required  in  the 
domain  D(a,b)  consisting  of  the  segment  a  <  x  <  to  of  the 
axis  y  =  0  and  of  the  part  of  the  half  plane  y  <  0  hounded  toy  two 
characteristics  of  equation  (Al)  passing  through  the  points  (a,0) 
and  (b,0) . 

t 

Chaplygin’s  equation  for  the  stream  function  as  a  function  of  the 
holograph  variables  can  he  brought  into  the  form  of  equation  (Al),  the 
domain  y<  0  corresponding  to  the  supersonic  region.  In  the  following 
considerations,  however,  no  use  will  be  made  of  the  special  form  of  the 
function  K. 

The  classical  theory  fails  to  treat  the  Cauchy  problem, 
equations  (Al)  to  (A3),  except  for  analytic  t  and  V,  because  the 
initial  data  are  given  along  the  .line  y  =  0  where  the  equation  ceases 
to  be  of  hyperbolic  type.  For  one  special  case,  that  of  the  Darboux— 
Tricomi’  equation 

+  ^yy  =  0  W 

the  solution  has  toeen  known  for  decades.  In  a  recent  paper, 
unavailable  in  this  country,  Christianovitch  treated  the  case  of  the 
Chaplygin  equation.  According  to  Frankl  (reference  1)  he  had  to 
assume  t(x)  and  v(x)  to  be  analytic  and  had  to  use  all  their 
derivatives  in  order  to  obtain  the  solution.  Frankl  (reference  l) 
considered  the  case  of  nonanalytical  Cauchy  data  'for  the  equation 

yf^  +  tyy  +  a(x,y)\|rx  +  b(x,y)\|ry  +  c(x,y)\|r  =  0  (A5) 

assuming  the  coefficients  a,  b,  and  c  to  be  analytic.  This 
treatment  is  based  on  special  properties  of  the  Riemann  function  for 
the  Darboux— Tricomi  equation  and  of  the  hypergeometric  function  and 
does  not  seem  to  lend  itself  easily  to  numerical  computations.  Note 
that  if  K (y)  is  analytic  and  K!(0)  f  0,  equation  (Al)  can  be 
brought  to  the  form  of  equation  (A5)  by  introducing 


as  a  new  independent  variable. 
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The  present  treatment,  while  limited  to  equations  of  the  form 
of  equation  (Al),  is  of  an  entirely  elementary  character,  requires  only 
very  simple  regularity  conditions  for  the  function  K,  and  remains 
valid  even  where  a  reduction  to  the  form  of  equation  (A5)  is  impossible 

(e.g.,  for  K  =  yl/3).  Besides  establishing  the  existence,  uniqueness, 
continuity,  and  differentiability  of  the  solution,  it  will  be  shown 
that  the  solution  depends  continuously  on  the  Cauchy  data  (as  noticed 
byFrankl),  as  well  as  on  the  coefficient  K.  These  continuity 
properties  imply  two  effective  methods  for  the  approximate  computation 
'of  the  solution:  Either  by  representing  it  as  a  series  of  particular 
solutions  of  the  Chaplygin  type  or  by  replacing  the  coefficient  K  by 
a  piecewise  constant  function,  in  which  case  the  determination  of  the 
solution  becomes  quite  trivial.  In  fact,  the  approximation  of  K(y) 
by  step  functions  constitutes  the  main  tool  of  this  investigation. 


2.  Statement  of  the  Problem 

The  function  K(y)  will  be  assumed  to  satisfy  the  following 
conditions : 

(1)  K( y)  is  defined  and  nondecreasing  for  y<  0 

(2)  K(y)  is  continuous  at  y  =  0 

(3)  K( y)  is  negative,  piecewise  continuous,  and  possesses  a 
piecewise  continuous  derivative  for  y  <  0 

Since  the  only  interesting  case  is  K(0)  =  0,  the  condition 
that  K(  y)  be  nondecreasing  involves  no  serious  loss  of  generality. 
(It  is  satisfied  in  the  case  of  Chaplygin's  equation.) 

A  function  t(x,y)  will  be  called  a  regular  solution  of  equa¬ 
tion  (Al)  if  \|r  possesses  continuous  partial  derivatives  of  the  first 
order  throughout  its  domain  of  definition,  and  piecewise  continuous 
partial  derivatives  of  the  second  order  satisfying  equation  (Al) . 

The  Cauchy  problem  mentioned  in  the  Introduction  may  be  stated  as 
follows : 

Problem  A.  Given  two  twice  continuously  differentiable 

functions  *r(x)  and  v(x)  for  a  <  x  ^  b.  To  determine  a 

regular  solution  ty(x>y)  of  equation  (Al)  defined  in  D(a,b) 

and  satisfying  the  initial  conditions,  equations  (A3) • 

It  will  be  convenient  to  consider  instead  a  more  general  problem 
(suggested  by  gas  dynamical  applications).  A  regular  solution  t  of 
equation  (Al)  determines  (except  for  an  additive  constant)  a 
function  cp(x,y)  such  that 
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<PZ  =  fy 
<Py  =  -&(j)Vx 


(A7) 


(Note  that  the  derivative  <py  might  possess  discontinuities.)  The 
complex— valued  function 


v(z)  =  <p  +  it 

► 


(A8) 


z  =  x  +  ly  j 


is  called  a  E-monogenic  function,  where 


£  =  Z(K) 


1  1 

1  K(y) 


(A9) 


Equations  (A7)  imply 


(A10) 


for  every  closed  rectifiable  curve  T  contractible  to  a  point.  It  will 
be  convenient  to  call  every  continuous  complex— valued  function, 
equations  (A8),  satisfying  equations  (A10)  a  E-*nonogenic  function. 

Part  of  the  theory  of  E-monogenic  functions  (references  12  and  13) 

( originally  developed  primarily  in  the  elliptic  case  K  >  0)  remains 


denotes  the  coefficientr-matrix  of  equations  (A7)  • 
that 


dx  -  K*  dy  =  0 


> 


L 


t  dx  +  q>  dy  =  0 
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valid  with  this  new  definition.  In  particular,  every  E-monogenic 
function  w  =  9  +  it  possesses  a  L-integral 

W(z)  =  f\(0  =  r  cp  dx  -  U  dy  +  i  M  dx  +  cp  dy  (All) 

J  Zq  «zo  t/zo 

W(z)  is  again  a  E-monogenic  function,  and 


(A12) 


Now,  if  is  a  solution  of  problem  A,  then  (for  a  proper 
normalization  of  cp) 


cp(x,0)  =  /  v(l)  d|,  a  <  x  <  b 

Ja 

This  remark  suggests  the  following  problem: 


(A13) 


Problem  B.  Given  two  continuous  functions  t(x)  and  v(x) 
for  a  <  x  <  b.  To  determine  a  E-monogenic  function  w  =  9  +  ity 
defined  in  D(a,b)  and  satisfying  the  initial  condition 

v(x)  =  J*  v(  I)  d|  +  1t(i),  a<  x<  b  (Al4) 


It  is  easy  to  see  that  if  w  is  the  solution  of  problem  B,  and 
if  i]r  possesses  continuous  first  and  piecewise  continuous  second 
derivatives,  then  is  a  solution  of* problem  A. 


3.  Statement  of  the  Results 

The  main  results  of  this  investigation  are: 

Theorem  1  (Existence  and  uniqueness  of  the  solution) . 

Let  r(x)  and  v(x)  be  arbitrarily  given  continuous  functions 
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defined  for  a  ^  x  <  b.  There  exists  one  and  only  one 
E(K)— monogenic  function  v  =  cp  +  ity  defined  in  D(a,b)  and 
satisfying  the  Cauchy  condition 


v(x)  =  J  v(!)  d|  +  1t(x),  a<  x<  b 


(A15) 


Theorem  2  (Continuous  dependence  on  the  initial  values). 

The  function  v  =  <p  +  i\|r  defined  in  theorem  1  satisfies  in  D(a,b) 
the  inequalities 


|cp|  <  AT  +  BN  I 


(A16) 


|t|  £  T  +  |yj  N 

where 

T  =  max  |i(x)J 
N  =  max  |v(x)| 


a  <  x  <  b 


A  =  A(y)  =  \/-K(y) 


B  =  B(x,y)  =  x  -  a  +  |y|A(y) 

Theorem  3  (Continuous  dependence  on  the  coefficient  of  the 
e quation) .  Let  Kn(y)  be  a  sequence  of  functions  satisfying 
conditions  (l)  to  (3).  Let  vn  =  cp^  +  ityn  be  the  Z(Kn)-Hnono genic 
function  satisfying  condition  (A15),  t  and  v  being  fixed 
continuous  functions.  If 
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lim  Kn(y)  =  K(y) 
n — 5^  00 


uniformly  for  0  >  y  >  -H,  where  H  is  the  least  upper  bound 
of  J  y  |  in  D(a,b),  then 


lim  wn(  z)  =  w(z) 
n — >°° 


uniformly  in  D(a,b),  w  being  the  E(K)-monogenic  function 
satisfying  condition  (A15).  (Cf.  the  remark  made  after  lemma  4.) 

Assuming  that  theorems  1  to  3  hold,  it  is  easy  to  establish 
similar  propositions  for  the  derivatives  of  the  function  w. 

•  Theorem  4.  If  r(x)  and  v(x)  are  continuously  differen¬ 
tiable,  then,  using  the  notation  of  theorems  1  and  2,  the  partial 
derivatives  cpx,  \(rx,  and  4y  exist  and  are  continuous  in  D(  a,b) 
whereas  <py  exists  and  is  continuous  in  D(a,b),  except  for  those 
values  of  y  for  which  K(y)  is  discontinuous.  Furthermore 

M  =  l^yl-  AT*  +  BN' 

\\\  <  T*  +  |y|  N* 

|cpy|  <  A3T»  +  A^Ii* 

where  T'  =  max  |t  '|  and  N'  =  max  |v  '  |  • 

Theorem  5.  If  r(x)  and  v(x)  possess  continuous 
derivatives  of  the  second  order,  then,  using  the  notation  of 
theorems  1,  2,  and  4,  the  partial  derivatives  cpxx,  txy> 
and  \yxx  exist  and  are  continuous  in  D(a,b),  whereas  the 
derivatives  cpxy,  cpyy,  and  tyy  exist  and  are  continuous 
in  d( a,b) ,  except  for  those  values  of  y  for  which  K(  y) 
or  K'(y)  is  discontinuous.  Furthermore 
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Krl  -  M  5  AT"  +  ™* 

I'M  ’  I VI  -a2t”  +  a2Ij|b" 

jcpy^  <  |E’(y)|  (AT*  +  |y|N')  +  a3t"  +  A2BN" 
where  T"  =  max  |t"|  and  N"  =  max  |v"|. 

■V 

Theorem  6.  Under  the  hypotheses  of  theorems  3  and  4  the 
partial  derivatives  of  the  first  order  of  9n  and  T]rn  converge 

(uniformly  in  D(a,b))  toward  those  of  9  and  .  Under  the 

hypotheses  of  theorems  3  and  5  the  partial  derivatives  of  the 

second  order  of  (fo  and  i|rn  (except  perhaps  q>n,yy)  converge 

(uniformly  in  D(a,b))  toward  those  of  9  and  The 

same  is  true  for  9n,yy  provided  that  K'n(y) ->  K'(y)  uniformly 

for  0  >  y  >  -H. 

The  proof  of  theorems  4  and  5  is  based  on  the  fact  that 
if  w  =  9  +  iijr  is  a  £-monogenic  function,  then  the  functions  $ 
and  if  defined  by 


®  +  i$  =  Jv  d^z 


possess  continuous  derivatives  (except  for  possible  discontinuities 
of  $y)  satisfying  equations  (A7),  so  that  if  the  functions  $ 

and  $  are  defined  by’ 

$  +  *1  =/(  $  +  if)  d£Z 

then  $  is  a  regular  solution  of  equation  (Al). 

Hence  theorems  4  and  5  follow  immediately  from  theorems  1  and  2  by 
applying  the  latter  theorems  to  the  Cauchy  problems  with  the  initial 
data  t  *  ( x) ,  v'(x)  and  r"(x),  v"(x).  The  proof  of  theorem  6 

proceeds  along  similar  lines  and  may  be  left  to  the  reader. 
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It  is  clear  that  the  functions  and  derivatives  whose  existence  is 
asserted  in  theorems  1/  2,  4,  and  5  exist  and  are  continuous  (or  piece¬ 
wise  continuous)  also  along  the  characteristics  hounding  D(a,b).  In 
fact,  the  functions  t(x)  and  v(x)  and  their  derivatives  have  been 
assumed  to  he  continuous  in  the  closed  interval  (a,b),  so  that  in 
order  to  establish  the  continuity  of  cp,  >|r,  and  so  forth  in  the 

closure  of  D(a,b)  it  is  necessary  only  to  continue  t  and  v  in  an 
appropriate  manner  over  the  interval  (a  —  e,  b  +  e) ,  €  >  0,  and 
apply  the  theorems  to  the  domain  D(a  —  €,  b  +  e).  Theorems  1  to  6, 
however,  remain  valid  if  T  and  V  and  their  derivatives  are  merely 
assumed  to  be  continuous  and  bounded  in  the  open  interval  (a,b). 

It  is  rather  obvious  what  modifications  would  be  necessary  if  T 
were  assumed  to  be  piecewise  continuous*  The  function  V  may  be 
assumed  to  be  bounded  and  measurable  without  endangering  the  validity 
of  theorems  1  to  3 7  as  will  be  seen  from  the  proof. 

4.  Formal  Powers 

A  significant- class  of  E-monogenic  functions  is  obtained  by 
setting  •  , 

Z(0)(z)  =1 

i’.z(°')(z)  =  i 


and,  for  n  =  1,  2,  .  .  .  , 

z(n)(z)  =  n  z(n~^)(£) 


i-z(n)(z)  =  n  J*  i*z(n  1)(C)  d£^ 


These  "formal  powers"  admit  the  representation  (see  reference  12) 
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z(”)(z)  ‘t.  (<”)l^iVvl(,) 
v=0 


(AIT) 


i*Z^n^(z)  =  i  (v)xn_vivY(V)(y) 

V  =0 


where 

Y(0)(y)  -  ?<0)(y)  -  1 


and,  for  n  =  1,  2,  .  .  .  , 


\ 

Y(2n+1)(jr)  =  (2n  +  J*  Y(2n)(n)  ^ 


Y(2n)(y)  =  2n  J* 

_ 

?(2n+1)(y)  =  ^  +  P 7  K^)~(2n)^)  dTJ 

Jo 


(Al8) 


^(Y)  =2nj[7Y(2n-1)(,)  dn  J 

Note  that  the  imaginary  parts  of  and  i*Z^n^  possess  continuous 

derivatives  even  if  K(y)  possesses  discontinuities. 

It  follows  from  the  preceding  formulas  that  (for  a  real  x) 
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V 


z(n)(i)  =  x11 
t.Z(n)(x)  =  lxn 


This  implies 

Lemma  1.  If 


M 

v(x)  =  21  15  ix^ 
3=0  J 


then,  for  every  real  number  a. 


aji-Z 


(3) 


(z) 


M  -h 

♦  Xr1 

3=0  J  + 


if 


( 3+1 ) 


(z)  - 


(A19) 


is  a  Z-monogenic  function  satisfying  the  initial  condition 


v(x) 


v(l)  d|  +  1t(x) 


\ 


Formula  (A19)  is  the  analog  of  D'Alembert’s  solution  of  the 
equation  of  the  vibrating  string.  In  fact,  if  K(y)  =  —A2  =  Constant, 
equation  ( Al)  becomes 


'*ix  2  ° 
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and  a  simple  computation  shows  that  in  this  case 

B<»)  =  g  {fX  +  +  ^X  "  +  2aEX  +  AY^  ”  ^X  “  Ay^nJ 

i*z(n)  =  ^J(x  +  Ay)n  -  (x  -  Ay)“]  +  |}(x  +  Ay)n  +  (x  -Ay)1*] 


so  that  the  imaginary  part  of  equation  (A19)  is  .  given  by 

*  =  ijx  ±  Ay)  +  r(x  -  AjQ  +  J_  P***7  v(|)  d| 
2  2Ac/x-Ay 


The  analog  of  Bernoulli's  solution  is  the  one  obtained  by  super¬ 
imposing  particular  solutions  of  the  Chaplygin  type.  (See  section  5  of 
the  ANALYSIS.) 

Consider  now  two  functions,  K^(y)  and  ^(y),  and  denote 

functions  (Al8)  formed  with  K  =  Kd  by  and  Yi^n^ 

where  i  =  1,  2.  Similarly,  let  the  E(Ki) -monogenic  formal  powers 

be  denoted  by  Zj^n)  and  i*Zj_(nt 

Lemma  2.  Set 


a  =  a(H)  =  max  |Kd(y)|  ,  0  >  y  -H,  i  =  1,  2 

&  =  5(H)  =  max  |%(y)  -  K^y))  , 

There  exist  constants  =  Mn'(a)  (depending  only  on  a 

and  n)  such  that 
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Yx(n)(y)  -  Y 2(n)(y)  I  <  5V  |yI n 

<  BM^*  |y|  n 

for  0  >  y  >  -H. 

The  proof  follows  easily  by  induction  on  n,  using  the  definition 
of  the  function  Y. 

Lemma  3»  There  exist  constants  Mn  =  Mn(a)  (depending  only 
on  a  and  n)  such  that  . 

|z1(n)(z)  -  ^(z)!  <  BMnlzl11 

1 

[i-Z^H^z)  -  i.Z2^n^(z)|  —  SMn  |z|  n 

for  0  >  y  >  -H.  Here  a  and  5  have  the  same  meaning  as  in 
lemma  2. 

The  proof  follows  immediately  from  lemma  2  and  from  equations  (A17) 


(n) 


(a) 


V“'(y)  -  V '(y) 


5*  Main  Lemma 

The  lemma  to  be  established  in  this  section  accomplishes  a  two¬ 
fold  purpose.  In  the  first  place,  it  shows  that  theorem  3  holds 
whenever  theorems  1  and  2  hold  for  the  functions  Kn  and  K.  In  the 
second  place,  it  shows  that,  if  theorems  1  and  2  hold  for  a  sequence  of 
functions  Kn  possessing  a  uniform  limit  K,  then  theorems  1  and  2 
hold  for  K.  Since  every  piecewise  continuous  function  is  a  uniform 
limit  of  piecewise  constant  functions  (step  functions),  it  follows 
that  theorems  1  to  3  are  valid  if  theorems  1  and  2  hold  for  the  case 
when  K  is  a  step  function. 

Lemma  4.  Hypotheses:  (a)  For  n  =  1,  2,  .  .  .  the  func¬ 
tion  "l^TyT”  satisfies  conditions  (l)  to  (3).  (b)  For  each  K^y) 
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theorems  1  and  2  hold,  as  well  as  similar  theorems  concerning  the 
Cauchy  problem  with  initial  data  on  the  line  y  =?  c  <  0. 

(c)  For  0  >  y  > -H  the  sequence  K^y)  converges  uniformly  to  a 
function  K(y)  satisfying  conditions  (l)  to  (3)* 

Conclusion:  Let  t(x)  and  v(x)  be  continuous  functions 
defined  for  a  <  x  <  b.  It  is  assumed  that  |y|  <  H  in  D(a,b). 
Let  wn(z)  be  a  z(Kn)~monogenic  function  satisfying  the  condition 


vn(x)  -  r  v(6)  dt  +  iT(x),  a  <  z  <  b  (A20) 

Then,  (a)  in  D(a,b)  the  sequence  wn  converges  uniformly  to  a 
(continuous)  function  w(z)j  (b)  w(z)  is  a  z(K)-monogenic; 

(c)  w(z)  satisfies  the  initial  condition  (A15);  (d)  w(z) 
satisfies  the  inequalities  (Al6);  and  (e)  w(z)  is  the  only 
function  in  D(a,b)  satisfying  assertions  (b)  and  (c). 

Remark.—  The  functions  wn  are  defined,  not  in  D(a,b),  but  in 
domain  Dn(a,b)  defined  in  the  following  proof.  In  view  of 
equation  (A22),  one  may  speak  of  the  convergence  of  the  sequence  wn 
in  D(a,b),  since  each  point  of  D(a,b)  belongs  to  all  but  a  finite 
number  of  the  regions  Dn(a,b) . 

Proof.—  At  first  it  will  be  assumed  that  assertion  (a)  holds. 

Set  wn  =  +  i>lrn  and  w  =  qp  +  iijr.  For  every  closed  curve  r 

in  D(a,b)  the  relations 


i 

s, 


<p  di  —  Kilf  dLy  =  0 


t  di  +  <p  cly  =  0 


(A21) 


must  hold  in  order  that  v(z)  Z(K) -monogenic.  Let  Dn(a,b)  denote 
the  set  defined  with  respect  to  Kn(y)  in  the  same  way  as  D(a,b)  was 
defined  with  respect  to  K(y).  Hypothesis  (c)  implies  that 


(A22) 


lim  Dn('a,b)  =  D(a,b) 

n— >  00 


ko 
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Hence,  for  sufficiently  large  values  of  n  the  curve  T  is 
in  Dn(a,b).  Then 


di  -  dy  =  0  1 


) 


dx  +  <pn  dy  =  0  J 


(A23) 


so  that  equations  (A2l)  follow  from  hypothesis  (c)  and  assertion  (a). 
Assertion  (c)  follows  immediately  from  equation  (A20).  Assertion  (d) 
is  true,  for  the  estimates,  expressions  (Al6),  hold  for  each 
function  wn  according  to  hypothesis  (a).  In  order  to  prove  the 

uniform  convergence  of  the  sequence  wn,  choose  a  positive  e  and 
determine  two  polynomials 

N 

t(i)  =  YI 

j=p 


M 

n(x)  = 

j=0  J 


such  that 


Jt(l)  -  t(x)|  <  6 
|n(x)  -  v(x)|  <  €  L 


(A2l+) 


a  <  x  <  h  J 

(This  is  possible  by  virtue  of  Weierstrass 1 s  approximation  theorem.) 
Let  Wn(z)  be  the  s(Kn)-monogenic  function  satisfying  the  condition 
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Wn(x)  =  J  n(l)  d|  +  it(x). 


kl 

(A25) 


Since  theorem  2  is  assumed  to  hold  for  K  =  Kn,  it  follows  from 
conditions  (A20),  (A24),  and  (A25)  that  for  all  points  in  Dn(a,b) 


|wn(z)  -Wn(z)|  <Le,  n  =  1,  2,  .  .  .  (A2 6) 

where  L  is  a  constant  independent  of  z,  n,  and  e .  Denote  hy 
and  i-Zn(m)  the  E(Kn)-monogenic  formal  powers  (formed  with 
respect  to  the  function  Kn) .  Since  Wn  is  uniquely  determined  hy 
condition  (A26) ,  it  follows  from  lemma  1  that 


Wn(z)  =  +  X  -  i-zJJ+l)(a)] 


Hence  lemma  3  and  the  uniform  convergence  of  the  sequence  Kn  imply 
the  existence  of  a  number  Mg  such  that  for  all  values  of  z  in  D(a,b) 

Kh-p(z)  -wm(z)|  <  e*  ®  >  Me,  p  =  1,  2,  .  .  .  (A27) 


By  virtue  of  expressions  (A26)  and  (A27) 

-  wm(z)|  <  (2L  +  l) e,  m  >Me,  p  -  1,  2,  .  .  . 

for  all  points  of  the  intersection  D(a,h)  DDm+p(a^) 

Since  L  is  fixed  and  e  arbitrary  and  since  Dn(a,b)  — ^  D(a,h) , 
assertion  (a)  follows. 

It  remains  to  prove  assertion  (e).  It  must  be  shown  that, 
if  g(z)  is  E(k) -monogenic  and  g(x)  =  0,  a  <  x  <  b,  g  vanishes 
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identically  in  D(a,b).  It  is  no  loss  of  generality  to  assume 

that  g(z)  is  continuous  in  the  closure  of  D(a,b)  (since  it  would  be 

b  —  e) ,  €  >  0) .  Then 


(x  +  ic)eD(a,b)  (A28) 


Furthermore,  no  loss  of  generality  is  involved  in  assuming  that  the 
real  and  imaginary  parts  of  g(z)  are  continuously  differentiable  (for 
otherwise  g(z)  could  be  replaced  by  its  e(k)— integral) .  The  argument 
leading  to  the  proof  of  assertions  (a)  to  (d)  maybe  repeated  to 
establish  the  existence  of  a  E(K) -monogenic  function  Gc(z)  satisfying 
the  condition  ' 


sufficient  to  consider  g  in  D(a  +  €, 
11m  max  |g(x  +  ic)|  =  0, 


Gc(x  +  ic)  =  g(x  +  ic ) ,  (x  +  ic)€D(a/b)  (A2<?) 

and  such  that  for  x  +  iy  =  zeD(a,b)  and  y<  c 


Gc(z)  =  lim  Gc  n( z ) 
n->oo  9 


(A30) 


^c,n(z)  being  a  z(Kn) -monogenic  function  with 

Gc,n(x  +  ic)  -  SU  +  i°),  (1  +  ic)eD(a,b)  (A31) 

Let  z0  =  x0  +  iyQ  be  any  point  of  D(a,b)  with  yQ  <  c.  It  follows 

from  hypothesis  (b)  and  theorem  2  that  there  exists  a  number  L  such 
that 


hc,n(zo)|  ^  L  max  |g(x  +  ic)| 
so  that  by  equation  (A30) 

|gc(zo')I  £  L  |e(x  +  ic)| 
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and  by  equation  (A28) 

Gc(zo)  =  0 

To  show  that  g  =  0  it  will  suffice  to  show  that  Gc(z)  =  g(z)  for  all 
values  of  c.  Set,  c  being  fixed, 

f(z)  =  g(z)  -  Gc(z) 

and 

F(z)  =  $  +  if 

where  d  is  such  that  d  +  ic  lies  on  the  characteristic  passing 
through  (a,0).  The  function  F(z)  is  z(K)-monogenic  and 

F(x  +  ic)  =0,  (1  +  ic)£D(a,b)  (A32) 

The  function  $  possesses  continuous  derivatives  and  fXy  and 

an  at  least  piecewise  continuous  derivative  $yy  It  satisfies 
equation  (Al)  and,  in  view  of  equation  (A32) ,  the  initial  condition 

f(x,c)  =  ?y(x,c)  =  0,  (x,c)eD(a,b)  (A33) 

But  for  y  <  c,  equation  (Al)  is  a  hyperbolic  type  and  the  classical 
theory  yields  the  result  $  =  0  in  D(a,b),  for  y<  c.  (The  fact 
that  K  and  K*  are  permitted  to  have  discontinuities  does  not 
impair  the  validity  of  this  argument.)  Thus  0=0  and, 
since  f  =  d^F/dx^,  g  =  Gc.  This  completes  the  proof. 
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6.  The  Step- Function  Case 

In  this  section  K(y)  is  assumed  to  be  a  step  function  (cf. 
figs.  5  and  6;  the  latter  shows  the  characteristics  and  the 
domain  D(a,b)). 

Set 


K(y)  =  -A/2,  7^ _±>7^7V  '  (A34) 


where 


0  =  yo  >  71  >  y2  • 
0  <A1  </^  .  .  . 


(A35) 


Theorems  1  and  2  will  be  proved  under  this  assumption.  In  view  of  the 
results  of  the  preceding  section,  this  insures  the  validity  of 
■theorems  1  to  3  in  the  general  case. 

Lemma  5*  If  K  is  given  by  equation  (A3^-)  >  a  complex— valued 
continuous  function  v  =  cp  +  ii{r  is  Z(k)— monogenic  if  and  only  if 
it  is  continuous  and  admits  the  representation 


v  =  Av  If  v(x  +  V)  “  SvC1  “  Avy)]  + 


i  |fv(x  +  Avy)  +  g^x  -  Avy)],  yy+1  >  7  ^  7y  (A36 ) 

(Formulas  connecting  fv  with  fv+1  are  given  in  the  following 
discussion.) 

Proof.—  It  is  easy  to  see  that  the  condition  is  always  sufficient 
and  that  it  is  necessary  when  ijr  possesses  continuous  derivatives  of 
the  first  order  and  piecewise  continuous  derivatives  of  the  second 
order.  Hence,  if  w  is  S-monogenic,  the  function 
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f(z) 


d.zz 


45 


admits  a  representation  of  the  form  of  equation  ( A36) .  Since  w  =  fxx 
the  assertion  follows. 

Lemma  6.  If  K  is  given  by  equation  (A3^) ,  theorem  1  holds. 
Proof.—  Set 


w  =  cp  +  it 


9  =  Av[fv(x  +  -  Mx  _Avyj] 


>_  > 


t  -  fv(x  +  \7)  +  gv(x  -  ^ )>  yv-i  “  y  -  y 


(A37) 


where 


flU)  =  |T(U  +  ~  I  v(**)  d?l 


1  Ua 


«1(S)  -IT(I)  fj  v(s,)  d5' 


(*38) 
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and  for  v  =  1,  2,  .  .  . 


fv+i(i)  =  |(x  +^)fv|j  -  (AV+1  -Av)yv]  + 

2^  “ “  ^v+1  +  ^v] 

«V+1^  =  l^1  "4^)fvf  +  ^+1  +  Av)yv]  + 

i(- +  ^)SvP  +  (Av+1  “  Av)yv] 


(A39) 


It  is  easy  to  check  that  w  is  £-monogenic  in  D(a,b)  and  that 
condition  (A15)  is  satisfied.  The  uniqueness  of  this  solution 
follows  from  lemma  5* 

Lenma_7.  If  K  is  given  by  equation  (A34),  theorem  2  holds. 

Proof.—  Because  of  the  linearity  of  the  problem  and  because  of  the 
special  form  of  inequalities  (Al6)  it  will  suffice  to  prove  these 
inequalities  under  the  assumption  that  at  least  one  of  the  functions  t 
and-  v  vanishes  identically.  Assume  first  that  v  =  0.  Using  the 
representation  of  the  solution  given  by  equations  (A37)  to  (A39) ,  it 
is  seen  that 

max  (|fv+l|  'I^V+ll)  ~  max  (|fv|  *|  ^vj)*  V  =  1,  2,  .  .  . 


and  that 
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Hence 


I'v+ll  -  T/2 


K+i|  s  T/2 


and  "by  equations  (A37) 


|<p  J  ^  AT 

M  ^  T 


since  for  K  given  by  equation  (3*0,  A(y)  >  Av  for  y  ^  Yv—i 
Now  suppose  that  t  =  0.  In  this  case  the  solution  admits  the 
rep  resentation 


where 


V-l 

v(z)  =  wv(z)  +  ^  wj(z),  yV-i  -  y  -  yv 


(A40) 


Wv  =  %  +  1%, 

vj  = 


(A4i) 


are  E-monogenic  functions  determined  by  the  initial  conditions 
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$1(x,0)  =  f  v  ^ 

I-l(x,0)  ho 

a  <  x<  b 

$v+l(x^v)  =  Mx^v). 

$v(x>yv)  *  0 

av  <  x  <  by,  v  =  2,  .  .  . 

<PV(I^v)  =  0 

Vv(x,yv)  =  $v(x^v) 

&v  ^  x  <  b^ ,  V  =  1,  2,  •  •  • 

Here  av  or  bv  denotes  the  abscissa  of  the  point  of  intersection  of 

the  line  y  =  yv  with  the  left  or  right  characteristic,  respectively, 
bounding  D(a,b).  Note  that  when  K  is  given  by  equation  ( 3U)  the 
characteristics  of  equation  (Al)  are  polygonal  paths  having, 

for  yv_i  -  y  -  yy>  the  sl°Pes  Ay-1  and  -Av-1.  It  is  easy  to 

verify  that 
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al  —  a  = 

ay  ”  aV  =  Ay(y^_ ^  7 y)  }  V  -  2,  3,  •  •  • 


For  0  ^  J  -  ?2_y 


|U,y) 


v(  5)  d?  + 


1  rx~Aiy 

2  Ja 


2AL  Jx—AjJ 


v(i)  d| 


v(|)  d| 


Hence 

|\(x,y)|  ^  (1  -  a)N 


(A42) 


(a43) 


bo  that. 


in  particular 


|\(ai*yi)|  -  (ai  -  a)K 


(A 44) 


and 


|*i,x(x'*l)| 


(A45) 


Also 


|f1(l^)|  <  iyiN 

For  yv-1  >  y  >  ^  ,  V  =  2,  3,  .  •  • 


(A46) 
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4  (x.y)  =  4 

V+Iv  ’  V 


(<\,7V)  +  2  jQ 


x+Av+l(y_7v) 


4y  x(S,yv)  d|  + 


1 

2 


j>*\+i(r-yv) 

°v 


®v,x(^yv)  dl 


r^i(x’7)  L 


x+A 


V+l 


(y-y,,) 


x^v+l(y~yv) 


\,x^yv}  dl 


Hence 

K+l(^y)l  -  |*v(vyv)|  +  Cz"av)  ““  |®v,xO'yv)l  (a4t) 
so  that,  in  particular, 

IViK+i'Wl  -  |Mwv)|  +  K+i  -  av)  (A48) 

and 

|4V+i,x(x^v+i)l  -  "“K'^^yv)!  (Al^9) 

Also  • 

|jv+l(x^)|  -  jy  -  yv|  max  |4  v,x(^yv)|  (A5°) 

From  inequalities  (A^5)  and  (A49)  it  follovs  that 

K,*(I,7v)|  ~  V  ~'^9  2,  .  •  .  (A5I) 
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Using  this  and  inequalities  (A4^)  and  (AW),  it  is  seen  that 

KXvyv)l  *  (*v  -  a)N 

so  that  by  inequalities  (A43)  and  (A47) 

|$v(x,y)|  <  (x  -  a)H,  y^  >  y  >  yv,  v  =  1,  2,  .  .  .  (A52) 

Also,  by  inequalities  (A46) ,  (A50) ,  and  (A5l) 

|$v(x,y)|  1  jy  -  Jy-i  -  yv*  v  =  1,  2,  .  .  .  (A53) 

so  that  in  particular 

^ K,  —  ^  =  2,  .  .  .  (A51*-) 

Using  this  last  inequality  and  the  previous  result  on  the  case  v  =  0, 
it  follows  that 

l^jj  -  AJy^  —  y J  =  l,  2,  .  .  . 

”  •  (A55) 
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Collecting  inequalities  (A52),  (A53),  and  (A55) ,  it  is  clear  that 
|cp(x,y)|  <L  jCx  —  a)  +  |y|A(y)jN  =  B(x,y)N 

|t(x,y)|  |y|N 

This  completes  the  proof. 
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A(M  <  1) 


r(M  =  1) 


Figure  1.—  Continuously  curved  arc  T  contained  in  boundary  of  domain  A 


M  <  1 


M  >  1 

1 1 11 II 


Figure  2.—  Flow  through  Laval  nozzle  with  transition  from  subsonic  to 
supersonic  flow  occurring  along  a  line  extending  across  the  nozzle 


cycloids 


cycloids. 
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Figure  8.—  Flow  through  Laval  nozzle  with  localized  supersonic  regions. 


M  <  1 


Figure  9*—  Flow  past  a  curved  surface  with  a  local  supersonic  region. 
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